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ABSTRACT 

We revisit an analytical model to describe the halo-matter cross-power spectrum and the halo 
auto-power spectrum in the weakly nonlinear regime, by combining the perturbation theory 
(PT) for matter clustering, the local bias model, and the halo bias. Nonlinearities in the power 
spectra arise from the nonlinear clustering of matter as well as the nonlinear relation between 
the matter and halo density fields. By using the "renormalization" approach, we express the 
nonlinear power spectra by a sum of the two contributions: the nonlinear matter power spec- 
trum with the effective linear bias parameter, and the higher-order PT spectra having the halo 
bias parameters as the coefficients. The halo auto-power spectrum includes the residual shot 
noise contamination that needs to be treated as additional free parameter. The term(s) of the 
higher-order PT spectra and the residual shot noise cause a scale-dependent bias function 
relative to the nonlinear matter power spectrum in the weakly nonlinear regime. We show 
that the model predictions are in good agreement with the spectra measured from a suit of 
high-resolution A^-body simulations up to k 0.2 /z/Mpc at z = 0.35, for different halo mass 
bins. 
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Clustering statistics of galaxies such as the two-point correlation 
function of galaxies or the Fourier-transformed-counterpart power 
spectrum are powerful tools to constrain cosmology. In particular, 
the baryon acoustic oscillation (BAO) experiment available from 
a wide-area galaxy redshift survey is recognized as a robust probe 
of cosmological distances. There are various on-going and planned 
galaxy surveys aimed at achieving high-precision BAO measure- 
ments over a wider range of redshifts: the SDSS-III Baryon Oscil- 
lation Spectroscopic Survey (BOSS)Q Subaru Prime Focus Spec- 
trograph (PFS) Survey F] and the ESA Euclid satellite experiment 


The BAO scale is one particular length scale measured from 
the pattern of galaxy clustering. Much more significant signal-to- 
noise ratios are inherent in the broad-band shape and amplitude 
information of the galaxy power spectrum at BAO scales. How- 
ever, to reliably use the amplitude information, we need to resolve 
various systematic uncertainties in the weakly nonlinear regime: 
the nonlinear clustering effect and galaxy bias uncertainty. There 
are promising developments towards a more accurate modeling of 
the nonlinear clustering of matter based on /V-body simulations 
Sprin gel et aLlp005);|Angulo et al.|p008);[T akahashi eTaT](f2009|; 
Nishimichi et al.|(|2009|>;|Valageas & Nishimichi|(|201 1|> as well as 



subara ] (|2008a); [T aruya & Hiramatsu (2008); [Saito et al.| ^2008); 
Nishi michi et al.|P007| >" 



perturbation theory (PT) of structure formation Juszkiewicz 1 1 98 1{>; 



Vishn iac|(|1983^|Makino et al.|(|1992|l;|J ain & Bertschinger (1994); 
Jeong & Komatsu| (2006| >; |Crocce & Scoccimarro| ( |2006| >; |MaT] 



The galaxy bias uncertainty is a harder problem, because phys- 
ical processes involved in galaxy formation/evolution are highly 
nonlinear and still very challenging to model from the first princi- 
ples. Hence a practical approach often used assumes an empirical 
parametrization of galaxy bias; the peak bias model|Kaiser|(|1984[>; 
|Mo & Whitel l [T996t ; |Sheth & Tormen| ( [T999l ; |Schmidt et al.|j2012) 
and the local bias model assuming a "local" mapping relation be- 
tween galaxy and matter distributions at each spatial posi tion|Coles| 
|l993t ; [Fry"& Gaztanaga| ( [T993) ; [Scherrer & Weinberg] ( |1998| >. In 
the peak bias model, a galaxy or more precisely halo is assumed to 
form at or around the peak of the initial matter density field, where 
a typical scale of the peaks corresponds to scales of the halo that 
host galaxies at low redshifts (although one halo can contain sev- 
eral galaxies inside). Thus the distribution of halos are by nature 
biased relative to the underlying matter distribution Kaiser ( |1984| >, 
because only the density peaks can be places to have halos today, 
while the under-density regions or the density minima are very dif- 
ficult (or impossible) to form halos. Further, the long-wavelength 
perturbation mode in the initial density field causes modulations in 
the heights of the small-scale peaks, and in turn alter subsequent 
formation of halos at low redshifts. The amount of halo bias can 
depend on the long-wavelength modes as a result of mode coupling 
in the weakly nonlinear regime relevant for BAO scales. The halo 



bias or peak bias models have been studied in the literature Mo &| 
White|fl996);|Sheth & Tormen|([T999);|Desjacques et al.| |207T] i; 
Scoccimarro et al.|(|2012^;|Schmidt et aL|<|2012), 



http://www.sdss3.org/surveys/boss.php 
~ http://sumire.ipmu.jp/en/2652 



3 http://sci.esa.int/science-e/www/area/index. cfm71areaid= 102 



The nonlinear effect on the galaxy power spectrum arises from 
two parts: the nonlinear clustering of matter (mostly dark matter) 
and the nonlinear bias relation between the galaxy and dark mat- 
ter distributions. At BAO scales in the weakly nonlinear regime, 
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we expect that the nonlinear clustering of galaxies can be accu- 
rately modeled by incorporating the PT of structure formation, the 
local bias model and/or the peak-background split bias model (halo 
bias). Such an attempt was first made i n|Heavens et al.|(|1998}, and 
followed by various works McDonald (20061; Matsubara (2008a); 



|Jeong & Komatsulp009l>; 



Saito et al.| (|2009);|M cDonald & Roy 



2009) ; |Manera et al.| ( |2010| >; |Baldauf et al.| ( |2010a| >; |Pollack et a" 
2Tm) ; |Sato~& Matsubara|pOTT) ; piareT7T1 ( (2rJT2l ); |Baldauf et al 
20121, which continuously show an improved understanding of the 



nonlinear galaxy power spectrum. 

In this paper, we revisit the method of modeling the nonlinear 
power spectra of halos, more explicitly halo-matter and halo-halo 
power spectra, by incorporating the PT, the local bias model and the 
halo bias. In doing this, we employ "renormalization approach" de- 
veloped in McDonald (2006) McDonald ( 2006 1 in order to re-sum 
contributions of the nonlinear matter clustering up to the higher- 
order loop corrections. This yields the term expressed by the prod- 
uct of the "full" nonlinear matter power spectrum and the renormal- 
ized linear bias parameter. The remaining terms, for which we keep 
the one-loop correction order based on the standard PT, give the ef- 
fect of scale-dependent halo bias in the halo power spectrum. Thus 
our approach is fully utilizing the recent improvement in modelling 
the nonlinear matter power spectrum (in this paper we will use the 
refined PT prediction developed in |Taruya~e t al. ( 2009)). We test 
the accuracy of our model predictions by comparing with the halo 
spectra measured from high-resolution A'-body simulations done in 
|Nishimichi & TamyalpOTT) . 

This paper is organized as follows. In Sec. [2] we develop a 
method of modeling the nonlinear halo-matter and halo-halo power 
spectra by incorporating the PT, the local bias model and the halo 
bias. In this section, we also show the detailed comparison of the 
model predictions with the simulation results for the halo catalogs 
of different mass bins. Throughout this paper we assume the A 
dominated, cold dark matter model (ACDM) as for our fiducial 
model which is consistent with Kom atsu et al.| ( |2011fr : the den- 
sity parameters of matter, baryon and the cosmological constant 
are Cl m0 = 0.279, Cl b0 /n m0 = 0.165, and Cl A = 0.721 (i.e. flat ge- 
ometry), the Hubble parameter h = 0.701, the tilt of the primordial 
power spectrum n s = 0.96 and the power spectrum normalization 
o- 8 = 0.817. 



2 PRELIMINARIES: PERTURBATION THEORY AND 
HALO BIAS 

Our model is based on three ingredients; the perturbation theory 
(PT) of structure formation Juszkiewicz] \ 198 1) ; Vishniac ( 1983| >; 



Fry] (19841); [Goroff et al.| (T986) ; |Suto~& Sasaki| fl99l]>; |Jain & 
Bertschinger (1994), the local bias model|Coles|(|1993^;|Fry & Gaz 



tanaga]i [7993[ ); |Scherrer & W einberg ( 1998 1 and the halo bias model 
Mo & White](|1996|>;|Sheth & Tormen|J1999[ l (also see |Cooray &] 
Sheth|2002||Bernardeau et al.|2002[ for thorough reviews). In this 
section, we briefly review the PT and the halo bias we will use for 
the following sections. 



2.1 Standard Perturbation Theory 

Throughout this paper, we consider a pressure-less, irrotational 
fluid system and assume cold dark matter as the dominant fluid 
component to drive gravitational instability of structure formation. 
The nonlinear dynamics in an expanding universe is fully charac- 
terized by the density fluctuation field, S,„, and the peculiar velocity 



field 9,„ Bernard eau et al.| p002). Given the initial conditions, the 
time evolutions of the fields are governed by the continuity equa- 
tion, the equation of motion and the Poisson equation. 

By using the standard PT, we can solve the nonlinear dynam- 
ics. The density fluctuation field at a given redshift z is expanded 



S m (k, z) 



d m (i)(k, z) + d m(2 )(k, z) + 6, n(3) (k, z) + ■ 



(1) 



The PT solution for the n-th order density fluctuation field can be 
found to be 



S m(n) (k, z) 



= Dl(z) J d\ 



1i ■■■d 3 q n 6 m([) (q ] )---S m([) (q„) 



XF„(q u --- ,q n )Sl 



*-2> 



(2) 



where 6(i)„, is the linear density field today, D + (z) is the linear 
growth rate normalized as D + (z = 0) = 1, and d 3 D (k) is the Dirac 
delta function. The n-th order density fluctuation field has the am- 
plitude of the order 0[(S(\ )m )"]. The growth rate can be computed 
by solving Eq. (7) in Oguri & Takada (2011). The Fourier ker- 
nel F„(q l , ■ ■ ■ ,q n ) gives a coupling between the different Fourier 
modes due to nonlinear clustering, and we will use the expression 
in Eq. (10a) of |Jain & B ertschinger]( |1994| >. Note that, although the 
form of the Fourier kernel is exact only for an Einstein de-Sitter 
model with Sl m = 1, it was shown to be a good approximation of 
the exact solution for a ACDM model. 



2.2 Halo Mass Function and Halo Bias 

Dark matter halos that host galaxies and/or galaxy clusters are use- 
ful tracers of large-scale structure, and can be used to infer the un- 
derlying dark matter distribution. However, the halo and dark mat- 
ter distributions are not the same, which leaves an uncertainty, the 
so-called bias uncertainty. In this paper we employ the halo bias 
model d eveloped in|Mo & W hite ( 1996); |Shefh & Tormen| fl999l > 
(also see |Cooray & Sheth|2002P 

Let us start with defining the halo mass function n(M, z)dM, 
which gives the comoving number density of halos in the mass 
range [M, M + dM] and at redshift z. We employ the mass func- 
tion given in Sheth & Tormen ( 1999 1: 



n(M)dM 



PmO 

~M 

PmO 



f(y)dv 

A[\ + (av)-"] 



Javexp 



dv 



(3) 



where p m0 is the mean mass density today; v = [<5 c /Z) + (z)<r„,(M)] 2 ; 
6 C is the threshold over-density for spherical collapse model; 
cr,„(M) is the present-day rms fluctuations in the mass density top- 
hat smoothed over scale R = (3M/4np m0 ) [l3 . We will throughout 
this paper employ the coefficients a = 0.75 and p = 0.3, which 
are obtained by comparing the fitting formula to A'-body simula- 
tions. Note that the normalization coefficient A is determined so as 
to satisfy the normalization condition j Q dvf(v) = 1. 

The mass function above holds only in an ensemble average 
sense, i.e. the average of the halo distribution over a sufficiently 
large volume such as the horizon volume. In other words, the num- 
ber density of halos in a smaller volume is modulated according 
to fluctuations of the underlying matter distribution within the vol- 
ume, <5,„. Employing the local bias model that the halo distribution 
is locally related to the underlying matter distribution, we assume 
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as the Fourier-space top-hat filter of smoothing scale R (also see 
|Smith et"aL 2007 ; Par frey et al.|201l] for the explicit expressions). 
The smoothing scale R is somewhat an arbitrary choice, and we 
naively expect that the smoothing scale needs to be chosen from 
the scales relevant for formation of halos of a given mass M. We 
have found that, on large scales we are interested in, the model pre- 
dictions without the smoothing well reproduce the nonlinear, scale- 
dependent halo bias seen from simulations. In this paper, therefore, 
we do not consider the smoothing function. 

In Fourier space, Eq. |4| can be rewritten as 



■5h(*,z) 



d 3 qi 



xs,„{qi,z) 



*-5> 

V i 

■S m (q„,z) + e(k). 



(7) 



Here, to have more generality as discussed in McDonald (2006), 
we added the term e(£) which gives the noise field uncorrelated to 
the matter density field, i.e. {eS m ) = 0. Also notice that the terms of 
(S" n ) in Eq. Q contribute only to the mode of k = in the Fourier 
coefficient, so are not relevant to the following results. 



Figure 1. Upper panel: the halo bias parameters, bi(M) and b 2 (M), as a 
function of halo mass, computed from Eq. |5J. We consider ACDM model 
and redshift z = 0.35. The nonlinear bias parameter, b 2 , is negative for 
low mass halos with M i. 2.5 X 10 13 , while it becomes positive and rapidly 
increases for the more massive halos. Lower panel: the renormalized PT 
prediction for halo bias functions for halos with masses M= 10 n ,10 13 and 
10' 3 " 5 M o //i (from bottom to top curves), respectively. We used Eqs. |l3| 
and to compute these curves. Note that we set the effective bias param- 
eter to = 0.9b[(M) (0.9 times the linear halo bias parameter) as implied 
from the simulations (see Sec. |4}, and set the residual shot noise term to 
5N = for simplicity. 



3 RENORMALIZED PERTURBATION THEORY FOR 
NONLINEAR HALO POWER SPECTRA 



In this section, we model nonlinear cross-power spectrum of matter 
and halos and nonlinear auto-power spectrum of halos by combin- 

iH2006) ; |Saito et al.| 



ing the "renormalized" PT approach 



McDonald 



(2009); |Jeong & Komatsu1j20(j9) ; [SaiTo et al.| ( |2dl lj with the local 



bias model, the halo bias and the perturbation theory described in 
the preceding section. 



that the number density fluctuation field of halos, <5h, is given as 

S h (x) = V -b n (M) [S" m {x) - <£(*)>] . (4) 

n 

The terms given as {6",) are added in order to satisfy the condition 
(<5h) = 0. The coefficients b„ are the halo bias parameters and given 
in terms of the derivatives of the halo mass function. To be more 
explicit, the 1st and 2nd-order bias parameters that are relevant for 
the following results are given as 



b 2 (M) 
where 

Cl 



1 + Ci + E, 



17\ 



2 1 1 - — I fa +Ei) + c 2 + E 2 , 



(5) 



1 



6 C 
2p/5 c 
1 + (av)i' ' 



C2 



av 

Ez 
E] 



(av-3), 

1 +2p 

6 C 



+ 2ci 



(6) 



When computing the power spectrum of halos including the one- 
loop corrections based on the standard PT, we formally need up to 
the 3rd-order bias parameter bj, but we will show that the contri- 
bution of £>3 is absorbed in the renormalized linear-bias parameter 
£>i ff . The upper panel of Fig. [T] shows the halo bias parameters, b { 
and b 2 , as a function of halo mass. 

Exactly speaking, Eq. |4| holds only for the smoothed fields, 
St,(x;R) and 6„,(x;R), where S(x;R) is defined by the convolu- 
tion d(x;R) = d(x) ® W(x;R) and W is the filter function such 



3.1 Halo-matter cross-power spectra 

First, let's consider the matter power spectrum defined as 

(S m (k)6 m (k')) = (2n) 3 P m (k)S 3 D (k + k'). (8) 

Using Eq. {T|l, we can find that the power spectrum including up to 
the one-loop corrections are given as 



P m (k; z) = z) + P m(l3) (t, z) + P m(2 2)(k, z), 



(9) 



where P)^(k;z) is the linear power spectrum, and P m (i3) and P„^ 22) 
are the one-loop corrections arising from the ensemble averages of 
the higher-order matter density fluctuation fields; <5,„(i)f5 m (3)) and 
(<5»i(2)^m(2)) ! respectively. The one-loop corrections at a given red- 
shift z can be computed once the linear power spectrum at the red- 
shift is specified: 



m(13) 



k 3 P];,(k ;z) 
252(2; 



tu) r 

n) 2 Jo 



drP L m (kr-z) 



12 



-42/ + -r(r 2 - l) 3 (7r 2 + 2) In 



158 + 100r 2 
1 + r 



■ ^r drpLjkr > z) 

X J d^P\-„(k^\+r 2 -2r^z) 



(3r + 7^- 10r^ 2 ) 2 
(1 + r 2 - 2rfi) 2 



(10) 



Similarly, using the standard PT and halo bias prescription, we 
can compute the cross-power spectrum between matter and halos 
of mass M, which is the quantity that halo-shear cross-correlation 
can directly probe. By inserting Eq. l|TJ into Eq. |7J, we can find 
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the formal expression of the halo-matter cross-spectrum in a self- 
consistent manner by including up to the one-loop correction terms 
of 0(^ a) ): 



P hm (k;M,z) 



b, + 



h + -b 2 



+bi [P,„(i3)(k) + P m{ 22)(k)] 



+bj 



J {2nf '" 



(q)P L J\k-q\)F 2 (q,k-q), 



(11) 



where a 2 = §d 3 ql{2n) 3 P\ n (q) and we employed notational sim- 
plification in the halo bias parameters; b, = bj{M). Thus, a for- 
mal implementation of the SPT based halo-matter spectrum yields 
the divergence term, i.e tr 2 ~ J q 3 P\ n {q)d\nq — > <x>, for a CDM- 
type power spectrum. In practice, since halo formation involves a 
smoothing of the underlying matter distribution as we discussed 
below Eq. j4|, the divergence does not arise in the power spectrum 
we actually observe. Also note that the prefactor coefficient of the 
linear power spectrum P\ n is independent of k. 

From the first two terms of the r.h.s. of Eq. we might 
re-write the two terms as 



68 

b\ H i,H b 2 

2 \ 21 



P]- n (k) + b,[P nm) {k) + P m(22) (k)\ 



Ibi + Sb l ]P h m (k) + b l SP m (k) 
[h x +Sb l ]P^m + 0(6b l SP m ), 



(12) 



where we have defined the notations Sbi = <r 2 (b 3 + 6Sb 2 /2l)/2 and 
SP m = P m (n) + P m (22), which are the higher-order contributions to 
the linear bias and the linear matter power spectrum by the order of 
0(<SLy) than leading order in the PT formalism. Motivated by the 
equation above as well as the similar idea proposed by McDonald 
|McDonald| ( [2006) , we propose the "renormalized" power spectrum 



P hm {k\M,z) = 
+b 2 (M) 



bfP^(k;z) 



r £q_ 

J (2tt)3 



P\M)P\,S\k-q\)F 2 {q,k-q). 



(13) 



Note that the terms in Eq. \\2\ are replaced with the nonlinear mat- 
ter power spectrum, />^ L , multiplied with the "effective" or "renor- 
malized" linear bias parameter, b'f as in McDonald ^2006} 



Thus 

the renormalized term can include the nonlinear corrections of mat- 
ter clustering up to the higher orders. From the PT viewpoint, this 
is not self-consistent in a sense that the term bf P^ includes the 



higher-order contributions than the one-loop order. Since |McDon-| 
|ald| ([2006 1 does not have the expression for the halo-matter cross 
correlation power spectrum, we will describe the difference be- 
tween our model and McDonald ( 2006 1 model once we derive the 
halo-halo auto correlation power spectrum in |3.2| 

There are several nice features in our renormalization pre- 
scription: 

• At the limit of small k or the linear regime, P\vn{k) — * 
b fP\„(k), because P^(k) -> P^ik) at the limit. 

• The effective linear bias b'f is a parameter, and is not related to 
the linear halo bias b\(M) due to the renormalization. Observation- 
ally, however, it can be determined by the cross-power spectrum at 
small k, e.g. measured from the large-scale signal of galaxy-galaxy 
weak lensing Mandelb aum et all| l |2012} . 

• The scale-dependent bias, relative to the nonlinear mass 
power spectrum, arises from the term that depends on the second- 
order halo bias b 2 (M) and the linear power spectrum and therefore 



is predictable once the background cosmological model, the halo 
masses and the redshift are given. 

We will test an accuracy of the renormalized cross-power spectrum 
(Eq. [T3| > by comparing the predictions with simulation results for 
halos of various mass ranges. We will below show how the use of 
the "full" nonlinear matter power spectrum in Eq. \\2>) can give a 
better fit to the simulation. 



3.2 Halo auto-power spectrum 

Similarly, we propose the renormalized auto-power spectrum of ha- 
los with masses M and M' (for generality of discussion, we con- 
sider the case that halos are in different mass ranges): 



P w (.k;M,M',z) 



bib\ 



o~ 
~2 



68 , 
b,b 3 + ^b x b 2 



68 

b[b' 3 + ^ h ib' 2 



P L m (k) + bib\ [P m(m (k) + P m(22) (k)] 



+ \ b ^S^ {q)p ^ k - q]) 

■^Lp\ n (q)P\„(\k - q\)F 2 {q, k-q) + SN 



bfbf'^ 1 

+ -b 2 b 
2 



P L Jq)Pl(\k-q\)-P L Jq) 2 



^4^(9)^(1* - ?l)^(9, k-q) + 6N\ 

(14) 

where we have introduced the effective linear bias bf and bf' for 
halos of masses M and M', respectively, and used the collapsed no- 
tations such as b 2 = b 2 (M) and b' 2 = b 2 (M') and similarly those for 
b[ and b' v The last term SN' denotes the residual shot noise term 
arising from (e(k; M)e*{k'\ M')) in Eq. jTJ as follows. As in the lit- 
erature, we refer to the term as the residual noise component after 
subtracting the Poisson shot noise, l/n h , without addressing the ori- 
gin. Assuming that the residual shot noise is constant over the scale, 
but may vary with halo masses, e(k, M) — > e(M), we will below 
study whether the model including the residual shot noise can give 
a better fit to the A'-body simulation result depending on the halo 
masses. In the second line on the r.h.s. of the above equation, we re- 
defined the shot noise term as SN' = SN + b 2 b' 2 J(d 3 q/(2nf)P] n (q) 2 
so that the term, which contributes to the scale-dependent bias, be- 
comes finite for the limit k — > <x>: 



2 j ^P L M)P\n(\k-q\) 

^b'tf ^[P L Jq)P L J\k- 



q\)- 



p^niqf] 



(15) 



Since Jd 3 qP^(q) 2 is constant, the residual shot noise term modified 
in this way is still constant in space, but varies with halo masses 
both through the dependence on b 2 (M)b 2 (M') and e. The constant 
SN' is treated as an additional free parameter for predicting Phh(k). 
The effective bias parameter bf in Eq. ( m4| is the same to that in 
Eq. d X 3 b up to the order Oid 2 ,^). Again we would like to stress 
some features of the renormalized halo power spectrum: 

• At the limit of small k or the linear regime, Pw(k) — > 
bfbf'P\ n (k) + SN'. 

• The scale-dependent bias depends on the linear power spec- 
trum and the halo bias parameters, b t (M) and b 2 (M), and therefore 
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is again predictable once the background cosmological model, the 
halo masses and the redshift are given. 

• The residual shot noise term SN needs to be included and 
treated as a free parameter. 

In McDonald (20061, b\ appeared in the third term of Eq. ( |14| is 
replaced with the renormalized Vf since b\ and Vf are approxi- 
mately same within the order of 0(S 2 m{l) ) while in our formalism, we 
keep the coefficients as halo bias, b\(M) as naturally derived from 
the PT and halo bias. We will below test the validity of Eq. i | 14} 
using simulations. In the following, we refer the residual shot noise 
term as SN, instead of SN', for simplicity. 

The lower panel of Fig. [T] shows the halo model predictions 
for the effective bias functions, which are defined in terms of 
the halo and matter power spectra: b cross (k) = Ph m (k)/P m (k) or 
b wt0 (k) = Phh(fc)/^hm(*) (see Eqs.[l3]and[n]>. Note that we did not 
include the residual shot noise contribution in this plot (set SN = 0). 
Our model predicts a scale dependent bias in the weakly nonlin- 
ear regime, and the degree of the scale dependence changes with 
halo masses. The top panel of Fig. [T] shows that bi(M) is negative 
for low mass halos, goes to zero around M =s 2 x 10 13 M e /h and 
then becomes positive for more massive halos. Hence, Eqs. jT3| 
and l |14[ l tell that the nonlinear bias due to bi(M) suppresses the 
power spectrum amplitudes for low mass halos, while it enhances 
for high mass halos. Also the model shows that, unlike the linear 
theory prediction, the scale-dependent halo bias generally differs in 
the estimators, Ph„JP m and Phh/Phm- Hence r t 1 in the weakly 
nonlinear regime, where r is the correlation coefficient of halo bias, 
defined as r = P hh / y/Ph,„P m - 



3.3 HOD model: Relating halos to galaxies 

Although we have focused on the halo power spectra, halos are 
not a direct observable from a galaxy redshift survey and need to 
be inferred from the distribution of galaxies. A practically useful 
approach to relate galaxies to halos is using the halo occupation 
distribution (H OD)|Peacock & Smith|pO^ ; [SeTjak| ( (2000l >; |Scoc-| 
|cimarro et al.H200l} . In this paper, for simplicity we assume that 



we can implement the method developed in Reid & Spergel Reid 



& Spergel (20091 (also see Reid et al. 2009, 2010, Hikage et al. 



2012a|b| l for reconstructing the halo distribution from the measured 
galaxy distribution. In the halo catalog, each halo hosts only one 
galaxy. 

In this setting, there is only one galaxy per halo. The cross- 
power spectrum of matter and galaxies and the galaxy auto-power 
spectrum are given in terms of the halo spectra as 



Pgm(k, Z) 



n CJ J 



JM^N H0D (M)P hm (k;M,z), 

>IM^N SW (M) 
dn 



dM' ^-N HO b(M')P w (k; M, W, z), 
dM' 



(16) 



(17) 



where the spectra P hm and P^> are given by Eqs. < | 1 3 [ > and \14) , 
respectively. Here A^hod(^) is the halo occupation distribution 
(HOD), and h g is the mean number density of galaxies, defined 
as hg = fdM (dn/dM)NuoD(M). Exactly speaking, the equations 
above are valid only if each galaxy is at the center of halo. When 
the galaxies have an offset from the halo center, we need to in- 
clude a convolution of the average offset distribution of galaxies, 
p oS (k;M), with the HOD distribution (see |Hikage et al.| < f2012a|b) 



for details). However, for the real-space power spectrum we fo- 
cus on in this paper, the effect is negligible on scales of interest; 
p fi(k; M) « 1 at the scales of interest. For the redshift-space power 
spectrum, the off-centered galaxies cause a significant Fingers-of- 
God effect. The main focus of this paper is the scale-dependent 
galaxy bias due to nonlinearities of structure formation, so we fo- 
cus on the real-space power spectrum. 

Recently Hamaus et al. Hama us~et al.| ( |201 \\ proposed a 
method, more directly based on the halo model approach, to model 
the nonlinear power spectrum of halos. In this method, the non- 
linear halo-matter and halo-halo power spectra are given by a sum 
of the linear power spectrum, multiplied with linear bias parame- 
ter, and the 1-halo term. We will below compare the performance 
of our method and theirs by comparing the model predictions with 
our own high-resolution A'-body simulations. 



4 TESTING THE NONLINEAR HALO SPECTRA WITH 
TV-BODY SIMULATIONS 

4.1 A'-body simulations and the halo catalogs 

To test the accuracy of our method for modeling the halo-matter 
and halo-halo power spectra in the weakly nonlinear regime, we 
use A'-body simulations for ACDM model done in Nishir mchi &| 



Taruya||201T|. In brief, we adopted 1280 3 A'-body particles and the 
box size of volume 1.5 (Gpc//z) 3 , and used the simulation outputs 
at z = 0.35. We defined halos using the Friends-of-Friends (FoF) 
finder algorithm with linking length 0.2 times the mean particle 
separation. For each halo, we use the total mass of member A'-body 
particles as the halo mass, i.e. the FoF mass, and the center-of-mass 
positions of the particles as the halo position. Then we computed 
the halo power spectrum from the discrete distribution of halos in 
each simulation realization using the cloud-in-cells interpolation 
method and Fourier transformation. To reduce the statistical scat- 
ters, we use the mean spectra from 15 realizations. To explore the 
validity of our model for different ranges of halo masses, we di- 
vided the halos into different mass bins. Table^shows the param- 
eters of the halo catalogs; we consider the halo catalogs divided 
into 9 mass bins, with mean masses ranging from 2.96 x 10 12 to 
7.03 x 10 13 M a /h. Note that, for the following results, we use the 
shot-noise-subtracted halo spectra, where we subtracted the theory- 
expectation l/n h from the halo spectra measured from the simula- 
tions (hi, is the mean number density of halos in a given simulation). 

The halo bias and the halo mass function we use are given as 
a function of the virial mass, rather than the FoF mass. We use the 
conversion relation M vir = 0.88M FoF at z = 0.35 to estimate the 
average virial mass for each mass bin using the method in |Hu &| 
Kravtsov (2003 1, where we assume that the FoF halo mass is close 
to the enclosed mass M lm , White (20021 inside which the mean 
density is 180 times the mean mass density and also assume that 
each halo follows an Navarro-Frenk- White profile Navarro et al. 
( |1997[ l with concentration parameter c vil = 4. We should note that 
this mass conversion only slightly changes the model predictions 
for the halo spectra, by up to a few % in the amplitudes. 

To compute the power spectrum for halos in the finite mass 
range used for the halo catalogs in Table [T] we use the following 
HOD in on our model (Eqs.|16|and|17[>: 



Mtod(M) : 



otherwise. 



(18) 
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Halo catalog 


bin 1 


bin 2 


bin 3 


bin 4 


bin 5 


bin 6 


bin 7 


bin 8 


bin 9 


AW10 12 [M Q /h] 


1.55 


2.18 


3.11 


4.37 


6.22 


8.77 


12.5 


17.6 


24.9 


M max /10 12 [M /h] 


4.86 


8.95 


15.3 


23.3 


35.4 


59.3 


104.4 


182.5 


- 


M h /10 12 [Mo/h] 


2.60 


4.08 


6.21 


8.22 


12.9 


19.1 


28.2 


40.6 


61.7 


n h /10" 4 [h 3 Mpc~ 3 ] 


15.7 


12.6 


9.46 


6.87 


4.87 


3.47 


2.43 


1.64 


1.09 


nnPrnik = 0.1 ft/Mpc) 


7.35 


6.78 


5.90 


4.97 


4.16 


3.57 


3.06 


2.52 


2.12 


nhP*h(k = 0.2 ft/Mpc) 


2.71 


2.54 


2.22 


1.88 


1.58 


1.34 


1.14 


0.93 


0.76 


fef (from P hm /P,„) 


1.06 


1.13 


1.22 


1.32 


1.44 


1.60 


1.78 


1.99 


2.27 


C (from SPT) 


0.66 


0.58 


0.42 


0.27 


0.13 


0.00 


-0.14 


-0.21 


-0.29 


C (from CPT) 


0.74 


0.65 


0.48 


0.32 


0.17 


0.04 


-0.10 


-0.17 


-0.24 


fel(M A ) 


1.17 


1.25 


1.34 


1.44 


1.55 


1.69 


1.91 


2.10 


2.33 


S 2 (M A ) 


-0.47 


-0.45 


-0.42 


-0.36 


-0.27 


-0.12 


0.22 


0.60 


1.19 



Table 1. Summary of the catalogs of simulated halos, built from the iV-body simulation outputs at z = 0.35 for ACDM model (see text for details). We use the 9 
halo catalogs, named as "bin 1", "bin 9", which are different in their mass bins; M m j„ and M max are the minimum and maximum masses to define each mass 
bin. Mi, and n/, are the average halo virial masses and the mean number density of halos in each halo catalog, respectively. As an indicator of the shot noise 
contamination to the halo power spectrum, we give the values of fih?hh at k = 0.1 and 0.2 /i/Mpc, respectively, measured from the simulations. The quantity 
is the renormalized linear bias parameter, which is estimated by comparing the PT model prediction to the simulation result for the halo-matter spectrum, 
Phm(k), at large scales (k < 0.05/i/Mpc) for each halo catalog (see text for details). The quantity C is the best-fit parameter to characterize the residual shot 
noise, estimated by fitting the PT model to the simulation result for the halo power spectrum, P\±(k), in the weakly nonlinear regime (see text). The values 
5i(M/,) and bi(Mh) are the halo bias parameters averaged by the halo mass function over the halo mass range. The effective bis parameter b c ® differs from the 
halo bias b\(Mh) by about 10% for the halo masses we consider. 



4.2 Halo-matter cross-power spectrum 

First, we study the halo-matter cross-power spectrum and define 
the bias function, for convenience of the following discussion, as 

Ph,n(k) 

(19) 

where P h .„, is the cross-spectra between dark matter (W-body) par- 
ticles and halos of the i-th mass bin, and P,„ is the matter power 
spectrum computed from the original A-body simulations. The shot 
noise is negligible for the cross-spectrum P^,,,. 

The data points in Fig.[2|show the bias function b^° ss (k) mea- 
sured from the simulations. The spectra P^m and P„, in each simula- 
tion share the same large-scale structure, and therefore the scatters 
due to the sampling variance mostly cancel in the ratio, yielding 
relatively smoothly-varying data points over k. The figure clearly 
shows that the halo bias has greater amplitudes for more massive 
halos. At sufficiently large scales or small k, k i 0.08 /j/Mpc, the 
halo bias appears to be constant, implying that the linear bias model 
is valid at the large scales. On the other hand, at the larger k, the 
simulation results manifest a scale-dependent halo bias for all the 
halo mass bins. 

The dot-dashed curves show our model predictions computed 
using Eq. \\2>) . To compute the model predictions, we need to fix 
one free parameter, the renormalized linear bias Uf , for which we 
determined Vf by fitting the model-predicted b cmss (k) to the sim- 
ulation b cross (fc) in the linear regime, at k < 0.05 /i/Mpc. Table [T] 
shows that the best-fit differs from the linear halo bias by about 
10%. Besides this free parameter, we used the input ACDM model 
parameters and halo mass range to compute the model prediction. 
Once the parameter Uf is determined, our model can be in remark- 
ably good agreement with the simulation results, including the k- 
dependence and the halo mass-dependence. For the largest mass 
bin, our model shows a sizable disagreement, possibly due to an 
inaccuracy of the halo mass function used for the model calcula- 
tion or the break-down of perturbation theory to describe too strong 
nonlinear bias. 

At large k, the perturbation theory ceases to be accurate, and 
is indeed not accurate enough up to k ~ 0.2 ft/Mpc. There have 



3 




0.1 0.2 0.3 

k [h/Mpc] 



Figure 2. The bias function £™ s (<0, defined as fe" oss (/t) = P hm (k)IP m (k) 
for halos of different mass ranges given in Table [l] The symbols are the 
simulation results, which are computed from 15 simulation realizations (see 
text for details), and clearly show a scale-dependent bias at k * 0.1 /i/Mpc. 
The red dot-dashed curves show the PT predictions including up to the 
one-loop correction (Eqs. |13| and |16[ . We determined the free parameter 
of the model prediction, the effective linear bias parameter , by fitting 
the model prediction to the simulation result up to k < 0.05 /i/Mpc for each 
halo mass bin. The scale-dependent bias is from a combination of the halo 
bias (bo) and the nonlinear matter power spectrum. The solid curves show 
the model predictions when using the improved PT model prediction given 
in Nishimichi & Taruya 1 2011 1 for the nonlinear matter power spectrum 
P™. L instead of the standard PT. 
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Figure 3. The halo-matter cross-power spectrum PhmM for halos of mass 
bin "4" (Table plotted in the unit of kP\ ml (k)ll7T 1 so that the features 
in the weakly nonlinear regime including the BAO features become promi- 
nent. The symbols are the simulation result as in Fig. [2] while the error 
bars are the statistical uncertainties at each k bins estimated from the 15 re- 
alizations. The dotted-dashed curve denotes the model prediction (Eq. |13) 
obtained by using the standard PT to compute the nonlinear matter power 
spectrum P„,(k) including up to one-loop corrections. The solid curve is the 
result if using the improved PT prediction (CPT) for P,„(k), which shows 
an improved agreement with the simulation result up to the higher k. To see 
the effect of scale-dependence bias arising from b% , the short dashed curve 
shows fej P^^ik), i.e the nonlinear matter power spectrum (CPT) multi- 
plied by the effective linear bias parameter. The dotted curve is the linear 
theory prediction, b c ^P\ n (k). For comparison, we also show the model pre- 
diction recently proposed in Hamaus et al. (201 1 1, which models the halo- 
matter power spectrum fully based on the halo model: the linear theory plus 
the 1-halo term given as bfP^(k) + P lh (k). 



been many efforts to improve the PT prediction of nonlinear matter 
power spectrum by including the higher-order loop corrections, e.g. 
the renormalized PT (Crocce & Scoccimarro 2006, Taruya et al. 
|2009| also see references therein). The solid curve shows the re- 
sults if we use the improved PT prediction for the nonlinear mat- 



ter power spectrum that is taken from Nishimichi & Taruya 



( |20 1 1 [ > (also see |Taruya et al.|2009) (the closure theory; hereafter 
CPT), instead of the standard PT including up to the one-loop cor- 
rection. The improved P^ n L has smaller amplitudes in the weakly 
nonlinear regime tiO.l h/Mpc at the redshift z = 0.35 than the 
standard PT predicts, yielding a slightly stronger scale-dependence 
of b{k){= PhmlP^) as evident from Eq. JT3I. The improved P^ L 



shows a similar-level agreement with the simulation results. 

In order to see the accuracy of our model for Ph„,(k), we com- 
pare the simulation result for P hm with the different model predic- 
tions in Fig. [3] Here we considered the intermediate halo mass bin 
(bin 4), and we use the parameter Vf determined in Fig. [T] En- 
couragingly, the figure clearly shows that our model prediction well 
agrees with the simulation Ph m (k) up to k ^ 0.2 h/Mpc, if we use 
the improved PT model (CPT) for the nonlinear matter power spec- 
trum. If we use the standard PT theory instead, the agreement is 
only up to k ^ 0.12 /j/Mpc, and the standard PT breaks down at the 



Figure 4. As in Fig. [3] but for the halo auto-power spectrum P^(k). 
The symbols show the simulation result for the halo catalog of mass bin 
4, where the standard shot noise term 1/% is subtracted from the mea- 
sured power spectrum. The dot-dashed and solid curves are the model pre- 
dictions (Eq. |14) , where we used the standard PT and the improved PT 
(CPT) to compute the nonlinear matter power spectrum P^(k), respec- 
tively. To compute the model predictions, we used the same linear bias 
parameter b c ^ in Fig. 2 and determined the residual shot noise parame- 
ter 6N from the fitting of each model prediction to the simulation result up 
to k = 0.2 /i/Mpc. To illuminate the scale-dependent bias, the dashed curve 
sho ws th e model prediction ignoring the terms of halo bias parameters in 
Eq. JT4J; i.e. (bf f ) 2 P% FT (k) + SN. The dotted curve is the linear theory pre- 
diction, {bffp\. 



larger k, in the weakly nonlinear regime. Hence we conclude that 
the apparent agreement of the standard PT up to the higher k bins 
in Fig.[T|is due to a cancellation of inaccuracies in the two spectra 
Phm and P m , in the numerator and denominator of the bias function 
b CIOSS (k). Comparing the solid and short-dashed curves shows the 
effect of nonlinear scale-dependent bias that arises from the term 
proportional to b 2 (M) in Eq. {13}. The scale-dependent bias be- 
comes important at k i 0.1 /z/Mpc for this redshift z = 0.35, and 
our model can reproduce the simulation result. 

The dashed curve shows the model prediction recently pro- 
posed in Hamaus et al. ( 201 1 1, where the nonlinear power spectrum 
is computed based on the halo model in combination with the halo 
bias parameters. The figure shows that the agreement is not as good 
as our model prediction. 



4.3 Halo auto-power spectrum 

The model parameters for the halo power spectrum Phh(k) are Vf 
and the residual shot noise parameter SN. For convenience of our 
discussion, we use the following parametrization of the residual 
shot noise term relative to the standard shot noise term: 



SN,,, 



(20) 



Fig. [4] compares the PT predictions and the simulation re- 
sults for Phh(&), as in Fig. [3] To compute the PT predictions, we 
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Figure 5. The bias function denned in terms of the halo auto-spectrum 
as b ant °(k) = y/Pbh(k)/P m (k) for each halo mass as in Fig. [2] The linear 
theory predicts no scale dependence and that the bias amplitude is the same 
to that of large scale limit for Ptm(k)/P m (k) in Fig. [2] The symbols are the 
simulation results for different halo masses, where the standard shot noise 
l/«h is subtracted. To compute the model predictions, we need to fix the 
free parameters: we used the same bf to that in Fig. 2 and determined 
the residual shot noise parameter C by fitting the modelprediction to the 
simulation result up to k < 0.2 /i/Mpc. 



determined the free parameters b'f and C for each halo catalog 
as follows. For b'f , we used the same values as those used for 
the halo-matter cross-power spectra in Fig. [2] For C, we deter- 
mined the value by fitting the model prediction to the simulation 
result up to fc raax = 0.15 or 0.20 /i/Mpc for the SPT or the im- 
proved PT (CPT), respectively. The maximum wavenumber fc max 
is chosen because each of the PT models for nonlinear matter 
power spectrum is accurate enough up to the fc max -wavenumber to 
within about 3% accuracy Nishimichi & Taruya < [201 1| >. In doing 
this fitting, we accounted for the statistical uncertainties in esti- 
mating the power spectrum from the 15 simulation realizations; 
i.e., we used the weighting in each £-bin, given as (AP hh ) 2 °c 
l/(2mt 2 A/t)[7 > hh (/t,) + l/h h ] 2 (fe, is the central value of the i-th k- 
bin and Ak is the width). The best-fit residual shot noise parameter 
C is about 30% compared to the standard shot noise for this halo 
mass bin (bin 4) as given in Table [T] 

As can be found from Fig. [4] our model prediction is in 
nice agreement with the simulation result, apparently up to k a 
0.25 ft/Mpc, if we use the improved PT prediction (solid curve) Q 
The nice agreement is found by including the residual shot noise 
contribution, which can account for a part of the nonlinear bias ef- 
fect. The standard perturbation theory cannot achieve the similar 
level agreement, even if varying the residual shot noise parameter. 
The figure also shows other model predictions, and the comparison 



4 Note that, exactly speaking, the improve PT (CPT) ceases to be accurate 
at k > 0.2 /i/Mpc, so the apparent agreement at the scales is as a result of 
the residual shot noise contribution, which happens to match the simulation 
result. 



of different model predictions manifests the importance of nonlin- 
ear clustering effect and scale-dependent bias in the weakly non- 
linear regime. Combining the results in Figs. [5] and [4] implies that 
the halo bias parameters £>i(M) and b 2 (M), in combination with the 
nonlinear matter power spectrum and the residual shot noise, can 
well reproduce the halo spectra Ph,„(k) and Phh(k)- 

Now we study another bias function defined in terms of the 
halo power spectrum as 



(k) 



*W*) 



(21) 



This bias function is different from the bias function £> cross (fc) we 
studied in Fig. [lias can be explicitly found from Eqs. \13\ and 
< | 1 4[ > (also Eqs. [16| and \T7\ . Fig. [5] compares the PT predictions 
and the simulations results for the bias function of each halo cat- 
alog b1" lo (k), where the PT predictions are computed by using the 
best-fit power spectra P hm and P hh as estimated in Figs. [2] and [4] 
The best-fit value of the residual shot noise parameter, C, for each 
halo mass bin is given in Table [T] The amount of the residual shot 
noise varies with halo masses, ranging from a few % to 70% com- 
pared to the standard shot noise term, and changes from negative to 
positive values from less to more massive halos. The bias function 
b mto (k) shows a scale-dependence over the range of k we consider. 
The scale-dependence of b ill "° indeed differs from that of b cmss (k) in 
Fig. [2] as the PT model predicts. The model predictions can fairly 
well reproduce the simulation results up to k m 0.2/i/Mpc, but then 
show a larger disagreement at the larger k than in Fig. [2] due to 
inaccuracies in both the model predictions for P hm and P hh . 



4.4 Comparison of halo bias model with other models 

The main difference between our method and the previous works 
McDonald| (2006)1; |Jeong & Komatsu] p009) l; [Saito et al.|j2009) ; 
Baldauf et al.H2010a| >; |Saito et al.| < |2011| >, which are both based on 
PT, is whether or not to incorporate the renormalized bias approach 
and the halo bias into the PT for computing the nonlinear power 
spectra. In this section, we compare our model predictions with 
other models that have similar forms. 

The nonlinear power spectra in these models are expressed by 
the following general forms (see Eqs.[T3]and[T4]for our model): 



PUk) ■■ 



a]P^(k) + 2a 3 P b2 (k) + ajP b , 22 + a 4 , 



(22) 



where 



PbMQ 



—^P m (q)P m (\k - q\)F 2 (q, k - q), 
= 2/(2^ [ P m(<l)Pm(\k ~ 9l) - P m (qf] 



(23) 



In terms of these forms, we can categorize the different models as 



Our model: c?i 



bf, a 2 



b 2 (M), a 3 = bi(M)b 2 (M), a 4 



6N, where b[(M) and b 2 (M) are the halo bias parameters and b'f 
and SN are treated as free parameters. 

• Saito et al. 2009: This is based on a renormalized perturbation 
theory originally proposed in McDonald (2006). Here the coeffi- 
cients are set to a { = b\, a 2 = b 2 , a 3 = b\b 2 , and a 4 = SN, and the 
three (bi,b 2 , 6N) are treated as free parameters to be determined by 
the fitting. This method is also studied in Baldauf et al.lpOlOb) . 

• McDonald 2006: This is similar to the method "Saito et al. 
2009", but uses the nonlinear matter power spectra P^f for P m 's 
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Figure 6. Comparison of the different PT-based model predictions for the halo-matter cross-power spectrum (left panel) and the halo auto-power spectrum 
(right panel), for the three halo catalogs of different mass bins (the 1, 4 and 8 mass bins in Table[T}. For illustrative purpose, the power spectra are normalized 
by the nonlinear matter power spectrum without BAO wiggles. For comparison, the dotted curves are the linear theory prediction (the power spectra in the 
denominator and numerator are both the linear power spectra). The different models are expressed by the similar forms (Eq. |22) , and have a different number of 
free parameters (2, 3 and 4 parameters), as indicated by legends and explained below Eq. J22J. The data with error bars at each k-bin are the spectra estimated 
from 15 realizations each of which has a 1.5 (Gpc/ft) 3 volume, and the error bar is the 1<t scatter of the central value at each k-bm, which is estimated by 
dividing the scatters of 15 realizations by Vl5; i.e. lcr statistical scatter for a volume of 15 X 1.5 (Gpc/ft) 3 . The best-fit parameters for each model are obtained 
by comparing the model predictions to the simulation taking into account the statistical errors (Eq. |24) . The shaded region at each k bin is the errors expected 
for a BOSS-like survey of 3.4 (Gpc//i) 3 volume, obtained by assuming that the error bars scales with survey volume (V,) as cr(P) oc 1 / V^s- Our model and 
the other models show a similar-level agreement with the simulations, to within the lcr error bars of BOSS-like survey up to k ^ 0.2 /j/Mpc, for different halo 
mass bins. 



in P b2 and Pt, t22 , instead of the linear spectrum. This method is in- 
tended to include renormalization for bias parameters as well as for 
the nonlinear matter power spectrum. 

• 4 free parameters: This is a variant of our model. The coeffi- 
cients are set to ct\ = Vf , a 2 = b 2 , «3 = b\b 2 , as, = 6N, and all the 
4 parameters (b'f, b t , b 2 , SN) are treated as free parameters. 

Note that, to have a fair comparison, we use CPT to compute P^, L (k) 
in the first term of Eq. j22} for all the above models. Thus the dif- 
ferent models have different ranges of their variations in the power 
spectra as a function of k with varying free parameters for a given 
cosmological model. 

Fig.[6]shows the different model predictions, described above, 
for Ph;« and for different halo mass bins, compared to the simu- 
lation results. To determine the free parameters in each model, we 
minimize the following x 1 by comparing the model prediction to 
the simulation result: 

^2 _ *y | ^psim _ pmodel^ £-1 ^psim _ pmodel^ (24) 

where P(Jq) = [Phm(kj), P\^(kj)], the power spectra with super- 
scripts "sim" or "model" are the simulated or model power spec- 
tra, respectively, C is the covariance matrix of the power spectrum 
computed from 15 realizations of the simulated power spectra, and 
C~' is the inverse matrix. We use the power spectrum information 
up to £ raax = 0.2 /j/Mpc. We included correlation between P hl „ and 



Phh at the same fc-bin in the covariance matrix, but ignored cor- 
relations between different k bins for simplicity Our choice of 
i mal = 0.2 /i/Mpc is based on the fact that the CPT prediction for 
P^, L (k) is accurate to within a 3% level up to k = 0.2 h/Mpc com- 
pared to the simulated spectrum at z = 0.35 as shown in Nis himichi| 
|& Tarayajj2TjTT] >. 

Fig. [6] shows that the different models well reproduce the sim- 
ulated P hm or P hh to an equal-level accuracy up to k ^ 0.15 or 
0.2 /i/Mpc in some cases. Again note that our model has least free 
parameters among these models, because our model uses the halo 
bias parameters b t (M) and b 2 (M) to model the scale-dependent bias 
of the nonlinear halo power spectra (Eqs.|22|(, and therefore restricts 
a range of the model variations compared to other models. Never- 
theless, our model appears to be reasonably accurate compared to 
other models. 

For comparison, the shaded region around the curves in Fig. [6] 
shows lcr statistical errors of the power spectrum measurements 
expected for a survey with volume coverage of about 3.4 (Gpc//?) 3 , 



5 This assumption would not cause any serious systematic errors as the 
scales we consider are in the weakly nonlinear regime and the non-Gaussian 
errors, which cause correlations between different k bins, are not significant 
as studied in |Takahashi et al.|p009| . 
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Phh when employing fc max = 0. 15 or 0.2 /j/Mpc, respectively) minus 
the number of free parameters (either 2, 3 or 4). Again our model 
and the models "Saito et al. (2009)" or "McDonald (2006)" are 
in a similar-level accuracy, and the model "4 free parameters" is 
slightly better due to more free parameters. Note that the reduced^ 2 
value is smaller than unity, partly because the central values of the 
simulated spectra are taken from the simulations of 22.5 (Gpc//i) 3 
volume and therefore the central value have smaller scatters than 
expected from a BOSS-like survey. 
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b halo 5N 

Figure 7. Comparison of the best-fit coefficients {cc\, aj, a^) in the 
model nonlinear power spectra (Eq. |22) with the best-fit parameters of our 
model, obtained from the 9 halo catalogs. Here we consider the different 
models as in Fig.[6](also see below Eq. |22) . Note that b\(M) and b2(M) in 
the jc-axis of the upper-right or lower-left panels are the halo bias parameters 
for different halo mass bins. The best-fit parameters are obtained by using 
the power spectra P^,,, and Phh up to k = 0.2 /i/Mpc for a BOSS-like survey, 
as in Fig.|6]or Table|2] The panels show that our model is almost equivalent 
to the models "Saito et al. (2009)" and "McDonald (2006)", implying that 
the halo bias parameters are a good approximation for the forms of the 
nonlinear power spectra given by Eq. \12) . 

which roughly corresponds to the volume of a BOSS-like survey 
with redshift range 0.5 < z < 0.7 and area coverage 10, 000 square 
degrees. We estimated the error bars by scaling the scatters at each 
k bin from the 15 simulation realizations assuming that the scatters 
scale with a survey volume as <t(P) oc 1/ \fV~ s . It can be found 
that our model is accurate at least enough up to k = 0.20 /i/Mpc 
within the \a errors for a BOSS-like survey, for a different range 
of halos masses. The linear theory is not accurate at BAO scales, 
at k i 0.1 h/Mpc, although the PT based model also ceases to be 
accurate at k t 0.2 /z/Mpc. 

Fig.^Jand Table|2]give a more quantitative comparison of our 
model with other models, where we consider the BOSS-like survey 
to compute the^ 2 given the expected measurement accuracies as in 
Fig. [6] Fig. ^compares the best-fit coefficients (art, . . . , a^), which 
can be read as effective bias parameters (see below Eq. [22| , with 
the results of our model; the halo bias parameters (i 1 (M),i 2 (M)) 
and the best-fit parameters Vf and SN. The figure shows that the 
results for "Saito et al. (2009)" or "McDonald 2006" reproduce 
the similar results to our model. The model "4 free parameters" 
shows a sizable difference from our result, especially for 0^3 and 
SN, implying that even small changes in the coefficients give the 
similar nonlinear power spectra for P^j and P hh at the scales. This 
also means a strong degeneracy between 0-3 and SN parameters. 
Table|2] gives the reduced,^ 2 of the best-fit model power spectra for 
each model and for each halo mass bin. Here we consider £ max = 
0.15 or 0.2 /j/Mpc for the maximum wavenumber to use for the 
model fitting. The degrees of freedom are defined by the number 
data point of the simulated spectra (58 or 78 in total for /> h „, and 



5 SUMMARY AND DISCUSSION 

In this paper, we have studied a method of modeling the nonlin- 
ear halo power spectra, P^ m {k) and Phh(k), by combining the PT 
approach of structure formation, the local bias ansatz and the halo 
bias. The nonlinearities of halo power spectra, which are devia- 
tions form the linear theory prediction, arise from the two effects: 
the nonlinear matter clustering and the nonlinear relation between 
the matter and halo density fields. In deriving the nonlinear halo 
power spectra, we employed the renormalization approach McDon-| 
|aldH2006) to re-sum the higher-order terms so that the terms are re- 
placed with the nonlinear matter power spectrum, multiplied with 
the "renormalized" linear bias parameter. The remaining terms in 
the nonlinear halo spectra are given as a function of the linear mat- 
ter power spectrum and the halo bias, and we kept the terms at the 
one-loop correction order. As a result, the halo-matter cross-power 
spectrum (Eq. 1 13^ is given as a function of the nonlinear and linear 
matter power spectra, the halo bias (biiM)), and one free param- 
eter, the renormalized linear bias parameter (ff ), which needs to 
be determined in the linear regime. Similarly, the halo auto-power 
spectrum (Eq.|14} is given as a function of the nonlinear and linear 
matter power spectra, the halo bias functions (bi(M), b 2 (M)), and 
the two free parameters, b'f and the residual shot noise parameter 
SN. Thus our method utilizes the recent development in an accurate 
model of the nonlinear matter power spectrum based on the refined 
perturbation theory and/or A'-body simulations. In our model, the 
halo power spectra are specified by cosmological parameters, halo 
mass, redshift, and a less number of free parameters. 

We showed that our model predictions for P\u„{k) and Punik) 
are in nice agreement with the simulation results, up to k = 
0.2 /j/Mpc, at simulation output z = 0.35 (see Figs. [3] and |4j, if 
using the improved PT theory prediction for the nonlinear matter 
power spectrum in the model calculation. The linear power spec- 
trum breaks down at k ^ 0.1 /?/Mpc. Thus our model might allow 
a factor 2 gain in the maximum wavenumber k m . dx up to which to 
include the power spectrum information when constraining cosmo- 
logical parameters. In the sampling variance limited regime, the 
factor 2 gain is equivalent to a factor 8 larger volume, yielding 
greater statistical power of the power spectrum measurement. In 
addition, a wider coverage of wavenumbers in the power spectrum 
amplitudes gives a higher sensitivity to some of intriguing cosmo- 
logical parameters such as the total neutrino mass |Saitoetal.| (2008 
2009, 2011) and the running index of the primordial power spec- 
trum. Thus developing a sufficiently accurate model of the nonlin- 
ear halo power spectrum is very important in order for us to have 
improved cosmological constraints, yet without having any signifi- 
cant biases in the derived parameters. 

Our model naturally predicts that, in the weakly nonlinear 
regime, the halo power spectra show a scale-dependent bias relative 
to the nonlinear matter power spectrum (see Figs.|2]and|5j. The PT 
based model naturally predicts that the two bias functions, defined 
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Nonlinear halo power spectrum 1 1 



Model 


Din 1 


Kin 1 

Din 2 


Kin 1 

Din j 


bin 4 


Kin ^ 

Din j 


bin 6 


Kin 1 

Din / 


Kin Q 

Din o 


bin 9 


This work 


0.38(0.54) 


0.52(0.75) 


0.51(0.75) 


0.62(0.86) 


0.53(0.70) 


0.60(0.76) 


0.32(0.34) 


0.23(0.33) 


0.46(1.3) 


Saito et al.2009 


0.34(0.43) 


0.49(0.62) 


0.45(0.78) 


0.68(0.80) 


0.67(0.61) 


0.30(0.31) 


0.32(0.35) 


0.23(0.26) 


0.50(0.95) 


McDonald 2006 


0.30(0.32) 


0.51(0.51) 


0.50(0.49) 


0.71(0.74) 


0.63(0.85) 


0.50(0.77) 


0.39(0.43) 


0.22(0.26) 


0.50(0.82) 


4 free parameters 


0.19(0.21) 


0.27(0.28) 


0.19(0.21) 


0.20(0.20) 


0.23(0.21) 


0.30(0.31) 


0.25(0.26) 


0.22(0.26) 


0.24(0.36) 


linear theory 


5.1 (13.4) 


5.4(15.7) 


4.3 (15.4) 


4.4(15.6) 


3.6(18.9) 


4.3 (23.3) 


8.4 (35.8) 


12.6(49.0) 


19.8(71.0) 



Table 2. Comparison of the different models as in Fig. [6] but the numbers in each row- and column are the reduced x 2 - values (xl) f° r tne best-fit power 
spectra of each model. The different columns are for the halo catalogs of different mass bins (see Table^. To compute the reduced^ 2 values, we obtained 
the best-fit model up to /: max = 0.15 /i/Mpc by fitting the model prediction to the simulation spectra assuming the errors for a BOSS-like survey, while the 
value in parenthesis is the value for £ max = 0.20 . The degrees of freedom for the^" fitting is: 58 or 78 in total for P^,,, and Phh f° r ^max 

= 0.15or0.2/i/Mpc, 

respectively, minus the number of model parameters (either 2, 3 or 4). Our model has 2 free parameters, the models "Saito et al. (2009)" and "McDonald 
(2006)" have 3 parameters, and the model "4 free parameters" has 4 parameters. The linear theory breaks down, but the different models of the nonlinear 
power spectra are equally acceptable for a BOSS-like survey. 



as b cmss (k) = P hm (k)IP m {k) and b* ut °(k) = y/P hh (k)/P m (k), differ in 
the weak nonlinear regime. Further, by incorporating the halo oc- 
cupation distribution (HOD) model, we can predict the nonlinear 
power spectra of galaxy-matter and galaxy-galaxy in the weakly 
nonlinear regime. We hope that our model can give a better de- 
scription of the nonlinear galaxy power spectra, and then allows for 
improved cosmological constraints via the measured power spec- 
tra. Our model can be further refined by including the higher-order 
loop corrections to the nonlinear bias functions. 

We showed that our model using the halo bias can give a 
good fit to the simulation results for the halo spectra. This offers a 
promising synergy between imaging and spectroscopic galaxy sur- 
veys, because a cross-correlation of the spectroscopic galaxies with 
images of background galaxies, the so-called galaxy-galaxy weak 
lensing, can directly probe the mean mass of host halos and in turn 
constrain the halo bias. Here the halo mass is constrained from the 
small-scale weak lensing signals arising from the mass distribution 
within one halo, which is complementary to the large-scale infor- 
mation of galaxy clustering at k i 0.2 /i/Mpc we focus on in this 
paper. This synergy is available if the spectroscopic and imaging 
surveys see the same region of the sky. This is the case for upcom- 
ing surveys: the BOSS and the Subaru HSC Survey, the Subaru 
HSC and PFS surveys, the Euclid, the WFIRST and a combination 
of the LSST survey with spectroscopic surveys. 

In this paper, we have focused on the real-space power spec- 
tra of halos or galaxies. Actual observable for galaxy redshift sur- 
vey is the redshift-space power spectrum of galaxies, which is af- 
fected by the redshift-space distortion effect due to peculiar mo- 
tions of galaxies. Towards a more accurate modeling of the non- 
linear galaxy power spectrum in redshift space, we need to further 
include the nonlinear coupling between the redshift-space distor- 
tion effect and the nonlinear galaxy bias. There are encouraging 
developments in modeling the redshift-space matter power spec- 
trum in redshift-space, based on the refined perturbation theory and 
A'-body simulations (Matsubara 2008b; Taruya et al. 2009 



2010 



|Tang et al.|20TT| |Matsubara|2011| |Sato & Matsubar a|2011 



. The 

redshift-space distortion effect due to virial motions of galaxies 
within halos, the so-called Fingers-of-God (FoG) effect, is harder 
to model, but Hikage et al. (2012a b) recently developed an empir- 
ical method to model the FoG effect based on the halo model and 
proposed a method to remove the FoG contamination by combining 
with galaxy-galaxy weak lensing measurement. It seems straight- 
forward to incorporate these methods in the method developed in 
this paper, in order to include all the effects, nonlinear clustering, 
nonlinear bias, nonlinear redshift-space distortion and FoG effect. 



This is our future work and will be presented elsewhere. These ef- 
forts are very important in order to attain the full potential of future 
high-precision galaxy surveys as well as to obtain unbiased, robust 
cosmological constraints from the surveys. 
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